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Abstract 

We investigate the finite temperature behavior of nonperturbatively improved clover 
fermions on lattices with temporal extent Nt = 4 and 6. Unfortunately in the gauge 
coupling range, where the clover coefficient has been determined nonperturbatively, 
the finite temperature crossover/transition occurs at heavy pseudoscalar masses and 
large pseudoscalar to vector meson mass ratios. However, on an Nf = 6 lattice the 
thermal crossover for the improved fermions is much smoother than for unimproved 
Wilson fermions and no strange metastable behavior is observed. 



Simulations with Wilson fermions suffer 0(a) scaling violations due to the 
dimension five operator that Wilson introduced to give the unwanted fermion 
doublers masses of order the cutoff, 1/a. These scaling violations are much 
larger than those in the glue sector, which are 0(a 2 ), and they can be numeri- 
cally quite large, necessitating use of small lattice spacings at large simulation 
cost to get results that can be reliably extrapolated to the continuum limit. 

The 0(a) scaling violations can be reduced to 0(a 2 ) by introducing another 
dimension five operator into the fermion action, the so called clover term, 



$sw Csw ft 



where T^ix) is a lattice transcription of the field strength tensor F^ u (x), 
usually taken from four open plaquettes looking like a clover leaf, as proposed 
by Sheikholeslami and Wohiert [1]. For the reduction of scaling violations 
to work the clover coefficient, c sw , needs to be determined nonperturbatively 
as a function of the gauge coupling. The ALPHA collaboration developed a 
method to do so within the Schrodinger functional framework [2]. For quenched 
QCD the nonperturbative clover coefficient is now known for gauge coupling 
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6/(? 2 > 5.7, corresponding to lattice spacings a ^ 0.17 fm [2,3]. A substantial 
reduction of scaling violations in this region from 0(a) to 0(a 2 ) has been 
verified nicely in [3]. 

The clover coefficient has recently also been determined by the ALPHA collab- 
oration for full QCD with two flavors of dynamical fermions for gauge coupling 
(5 = 6/g 2 > 5.2, corresponding roughly to lattice spacings a < 0.14 fm [4]. To 
be precise, the clover coefficient was determined for (3 > 5.4 and fitted to a 
ratio of polynomials in g 2 . At j3 = 5.2 a numerical consistency check with the 
value extrapolated with this function of g 2 was performed (see [4] for details). 
Preliminary first results of hadron spectroscopy and the heavy quark poten- 
tial in dynamical simulations with this nonperturbatively determined clover 
coefficient have appeared in [5]. 

Arguably the largest lattice artifacts in simulations with dynamical Wilson 
fermions have been observed in simulations probing the finite temperature 
behavior in the vicinity of the deconfinement /chiral symmetry restoration 
transition or (at finite quark mass) crossover. The most likely scenario for 
the behavior of two-flavor QCD at finite temperature is a rapid crossover 
for finite quark mass, turning into a second order chiral symmetry restoring 
phase transition in the massless limit [6] . However, simulations with dynamical 
Wilson fermions showed strange, unexpected behavior, e.g. first order phase 
transition like signals at intermediate quark masses, that softened again at 
smaller quark masses [7,8]. 

It has been argued [9] that this strange and unexpected behavior is due to 
effects of the Wilson pure gauge action in its so called "crossover region", 
where the plaquette varies sharply with changes in the gauge coupling, which 
feeds back to the fermions, rather than being due to artifacts in the fermion 
action. Indeed, their simulations did not show any evidence for first order like 
signals. Similarly, a study using both an improved gauge and a clover improved 
fermion action, with so-called tadpole improved coefficients, found a smoother 
behavior in the thermal crossover region than the simulations with unimproved 
Wilson action for both gauge and fermion sectors [10]. From the point of 
view of the Symanzik improvement program both these simulations still have 
0(g n a) errors of unknown magnitude. Furthermore, it is not clear, in the 
study where both gauge and fermion action are improved, which improvement 
is more important in smoothening out the thermal crossover behavior. 

In this letter, we study the effect of the nonperturbative improvement of the 
Wilson fermion action on the behavior at the finite temperature crossover. 
We are interested in simulations with small lattice extent in the temporal 
direction, i.e. at large lattice spacing, where the simulations are relatively 
cheap. 
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Table 1 

The clover coefficients, c sw , used in the simulations and estimates of the k c 's. They 
were obtained from Ref. [4]. 
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Fig. 1. (ReP) for the simulations with iVj = 4 and volume 8 3 at /3 = 5.4 (squares) 
5.3 (diamonds) and 5.2 (octagons) and with volume 12 3 at f3 = 5.2 (crosses). 

The largest coupling for which the nonperturbative value of c sw is known is 
(3 = Q/g 2 = 5.2. We performed simulations on 8 3 x 4 lattices for (3 = 5.4, 5.3 
and 5.2 and various values of k in the thermal transition/crossover region. The 
values for c sw used, obtained from Ref. [4] are listed in Table 1. Measurements 
were taken, after thermalization, over 500 trajectories away from the crossover 
region, and over up to 4000 trajectories in the middle of the crossover region. 
We show in Figure 1 the real part of the Polyakov line expectation value and 
in Figure 2 the average space-like plaquette. A clear crossover is seen for all 
three gauge couplings, becoming sharper as the coupling is increased (as (3 is 
decreased). For the largest coupling, (3 = 5.2, we also simulated on a larger 
spatial volume, 12 3 x 4. As can be seen from the figures, there is no evidence 
for finite volume effects. 

The thermal crossover appears quite rapid, reminiscent of unimproved Wil- 
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Fig. 2. Space-like plaquettes, (TrU Pa /S), for the same lattices and with the same 
symbols as Fig. 1. 

son simulations [7], as compared to thermodynamics with a tadpole improved 
clover fermion and Symanzik improved gauge action [10]. But, of course, com- 
parisons should not be made in terms of bare parameters, such as k. We 
have therefore computed hadron masses and the heavy quark potential at the 
thermal crossover points. For /3 = 5.2 hadron masses were also computed on 
both sides of the thermal crossover. Most of these measurements were done 
on 8 3 x 16 lattices. In one case, for (3 = 5.2 and k = 0.1340, where the meson 
masses are lightest, we also have preliminary results from a simulation on a 
larger 16 3 x 32 lattice. There, we do see finite size effects on the masses of 
about 10%, but much less on the mass ratio. We suspect that the masses for 
the next lightest quark mass, at (5 — 5.2, k = 0.1330 are also affected by 
some small finite size effects, while for the other cases the finite size effects are 
expected to be smaller than the statistical errors. All the results are collected 
in Table 2. 

Even at the strongest coupling, for which the nonperturbative clover coef- 
ficient is known, the thermal crossover for N t = 4 lattices occurs at heavy 
pseudoscalar meson mass and large pseudoscalar to vector meson mass ratio. 
We therefore also considered thermodynamics on an N t — 6 lattice at (5 — 5.2. 
(ReP) and (TrU Ps /3) are shown in Figures 3 and 4, where they are compared to 
the results from N t = 4. The crossover seems somewhat smoother for N t = 6, 
and it occurs at lighter meson masses. However, at 0.85 the pseudoscalar to 
vector meson mass ratio is still rather large. 
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Fig. 3. (ReP) for the simulations at (3 = 5.2 for lattice size 8 3 x 4 (octagons), 12 3 x 4 
(crosses) and 12 3 x 6 (diamonds). 
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Fig. 4. Space-like plaquettes, (TW Ps /3), for the same lattices and with the same 
symbols as Fig. 3. 
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Table 2 

Masses of pseudoscalar and vector mesons and their ratio, and values for the string 
tension from L 3 x (2L) lattices. Points marked by an asterisk lie on the N t = 4 
crossover, and the point marked by a plus on the N t = 6 crossover. 
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Comparing Figure 1 with similar plots for unimproved Wilson fermions such 
as Fig. 5 of [7] the crossover for the improved clover fermions appears to be 
even more pronounced than for the unimproved Wilson fermions. Comparing 
as function of the bare hopping parameter k, though, can be misleading and 
we therefore follow the strategy of Ref. [10] and make a comparison as function 
of the quark mass, or equivalently as function of m 2 PS . From Ref. [11] we can 
see that the crossover for Wilson fermions at f3 — 5.12 and k = 0.1700 occurs 
at a pseudoscalar to vector meson mass ratio m PS /m v = 0.899(4). It appears 
that the unimproved Wilson fermion data at /3 — 5.1 of [7] are the set that 
most closely matches ours at (5 — 5.2. In addition, from Ref. [12] we have some 
mass measurements in the thermal crossover region. While the mps/my ratios 
at the thermal crossover are comparable, the masses in lattice units are quite 
different. We therefore decided to plot (ReP) as function of (mps/my(KT)) 2 , 
where my(Ky) is the vector meson mass at the thermal crossover point kt- The 
comparison is shown in Fig. 5. It appears that the crossover for the improved 
fermions with Nt = 4 is somewhat sharper than for the unimproved fermions. 

The crossover for the improved fermions on the N t = 6 lattice, on the other 
hand, is much smoother. The pseudoscalar to vector meson mass ratio at the 
crossover corresponds to that for unimproved Wilson fermion at (5 — 5.22, 
k = 0.17 [8]. This is the region where first order like metastable states were 
observed in [8]. So here the improvement seems to help. 

In conclusion, we have carried out finite temperature simulations with non- 
perturbatively improved Wilson fermion in the region of largest coupling, and 
hence largest lattice spacing, for which the nonperturbative value of the clover 
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Fig. 5. Comparison of (ReP) for the simulations with clover improved and 
non- improved Wilson fermions at similar mpsjmy ratios at the thermal crossover 
as function of (mps /my («t)) 2 , giving a comparable measure of the quark masses, 
for N t = 4. Also shown are the clover results for N t = 6. 

coefficient c sw is known. The thermal crossover on Nt — 4 lattices occurs at 
very heavy pseudoscalar meson masses and large pseudoscalar to vector meson 
mass ratios. The crossover appears somewhat sharper than for unimproved 
Wilson fermions at comparable mps/my ratios. However, for the improved 
fermions the masses (in lattice units) are considerably larger, and the thermal 
crossover could still be significantly influenced by the deconfinement transi- 
tion in the pure gauge theory. While the thermal crossover on an N t — 6 
lattice at the strongest accessible coupling still occurs at a large mps/m v 
ratio, it has become much smoother, in particular compared to unimproved 
Wilson fermions, which show strange first order like behavior for comparable 
mpsfmy ratios. To get the thermal crossover to occur for smaller mps/my 
ratio, one needs to use lattices with larger temporal extent N t or try to push 
the nonperturbative determination of c sw to stronger couplings. 
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